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NONINVASIVE NONLINEAR SYSTEMS AND METHODS 



FOR PREDICTING SEIZURE 



This invention claims priority to US Provisional Application Serial Number 
5 60/410,695 filed on 13 September 2002. The entire disclosure of the priority application is 
specifically incorporated herein by reference. 

This invention was supported in part with grant ROl NS0368O3 fi-om the National 
Institutes of Health. The United States government may have rights in this invention. 

10 FIELD OF THE INVENTION 

The present invention relates to methods and devices for noninvasive nonlinear 
prediction of ictal onset in patients afflicted by neurological disease. In particular, the 
present invention provides methods and devices for noninvasive nonlinear prediction of 
seizures in patients afflicted with epilepsy. The devices and niethods preferably being based 
15 on analysis of two or more electroencephalogram (EEG) recordings, one set of recordings 
taken fi'om an electrode close to the region of ictal onset, and a second or more set of 
recordings (e.g., concurrent readings) taken firom a region remote firom the region of ictal 
onset. 

20 BACKGROUND OF THE INVENTION 

Epilepsy afTects between 40 and 50 million people worldwide — ^nearly 1% of the 
world's population. In the U.S. epilepsy affects about 2.5 million individuals at any given 
time with about 150,000-200,000 new cases diagnosed per year. (Begley et al.. Epilepsia, 
41:342-351 [2000]). Due to the persistent stigma surrounding neurological diseases such as 

25 epilepsy, the actual number of people affected by the disorder may actually be even higher. 

In approximately 75% of affected individuals, the epilepsy has no identifiable cause. 
In the other 25% of afflicted individuals, the most common causes of epilepsy are brain 
injuries and trauma {e.g., hypoxia/anoxia, especially during birth), substance abuse, 
infections (e^., meningitis), brain tumors or stroke, genetic defects, and degenerative 

30 disorders such as Alzheimer's disease. 
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An q)ileptic incident is characterized by intermittent interruption of normal brain 
function by sudden and often intense periods of synchronous neural discharge, resulting in 
either convulsive seizures, or more subtle alterations in neurological function such as brief 
lapses in full consciousness. 

5 Epilepsies come in many forms and have many causes. A central distinction in the 

classification of epilepsies is between those with origins in an identifiable region of the 
brain, called focal or partial epilepsies (e.g., mesial temporal lobe epilepsy), and those with 
no well defined site of origin, called generalized epilepsies. Each of these two major 
categories may be subdivided into symptomatic (due to a known cause e.^., a brain tumor) 

1 0 type or idiopathic (of unknown cause but commonly suspected to be heredofamilial in 
origin) type. Symptomatic epilepsies are often medically intractable. A diagnosis of 
intractable epilepsy usually comes after two or three years of unsuccessful treatment with 
standard anti-seizure medications in a compliant patient. 

While it is useful for neurologists to try to group epilepsies into categories for 

15 treatment, in reality, there can be significant overlap in the physical manifestations of a 
seizure disorder between patients with different types of epilepsy. One thing is clear, 
different epilepsies result from malfunctions in particular areas of the brain. Medicine has 
yet to provide well defined parameters for the different types of epilepsies, and has chosen 
instead to adopted a landscape classification approach in which various features such as 

20 clinical history, interictal (between seizures) EEG manifestation and results of neuroimaging 
(e.^., MRI) serve as useful landmarks. To date there is essentially no understanding of why 
a particular patient will have a seizure at any point in time. As a consequence, the 
prediction and treatment of epileptic seizures in individual patients remains very 
challenging, 

25 The unpredictability of ictal (seizure) onset in an individual afflicted with epilepsy is 

perhaps the most difficult aspect of Uving with an epileptic condition. Many types of 
epileptic seizures totally incapacitate the afflicted individual for from moments to hours. 
Loss of one's facilities while controlling a vehicle or operating machinery, among many 
other things, can lead to potentially dangerous situations for the epileptic and for others at 

30 large. 
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What are needed are new portable and noninvasive methods and devices for reliably 
predicting ictal onset in a variety of patients. Such methods and devices will allow those 
suffering from many types of epilepsies heretofore unknown freedoms. 



5 SUMMARY OF THE INVENTION 

The present invention relates to methods and devices for noninvasive nonlinear 
prediction of ictal onset in patients afflicted by neurological disorders. In particular, the 
present invention provides methods and devices for noninvasive nonlinear prediction of 
seizures in patients afflicted with epilepsy. The devices and methods preferably being based 
10 on analysis of two or more electroencephalogram (EEG) recordings, one set of recordings 
taken from an electrode close to the region of ictal onset, and a second or more set of 
recordings (e.g., concurrent readings) taken from a region remote from the region of ictal 
onset. 

Embodiments of the present invention may be configured to predict the onset of ictal 
15 episodes in subjects having partial epilepsies including, but not limited to, benign occipital 
epilepsy (benign focal epilepsy with occipital paroxysms), benign rolandic epilepsy (benign 
focal epilepsy with centrotemporal spikes), frontal lobe epilepsy, occipital lobe epilepsy, 
mesial temporal lobe epilepsy, other forms of temporal lobe epilepsy, and parietal lobe 
epilepsy. 

20 Broadly speaking, there are two types of epilepsy, focal and generalized. In focal 

epilepsies there is thought to be a specific region of the brain from which the seizures 
originate (although, as we shall discuss below, this notion is vague and the reality may be 
considerably more complex). The most common type of focal epilepsy is temporal lobe 
epilepsy (TLE), in which the region of ictal (seizure) onset is in one (rarely both) of the 

25 temporal lobes. Generalized epilepsies are those in which there is no clearly identifiable site 
of ictal onset. 

Of all epilepsies, about 50% are focal epilepsies, and of these roughly 70% are 
epilepsies of the temporal lobe. Of patients with focal epilepsy, roughly 25% suffer a 
medically refractory condition, so that the only possible treatment currently available to 
30 them that might result in control of their seizures is surgical resection of part of the temporal 
lobe. For these patients, in particular, the present invention provides a reliable portable 
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method for seizure prediction. The present invention allows the patient to reliably position 
himself in a safe environment (e.g., not driving, away from machinery, etc.) to weather the 
seiziire. Additional embodiments of the present invention further incorporate one or more 
devices (e.g., electrical stimulation, medication dispensers, and the like) for administering 
5 therapies sufficient for aborting ictal onset or for lessening its effects. 

While the present invention is principally directed to providing methods and devices 
for predicting the onset of epileptic seizures arismg at a particular focal point in the brain, 
additional embodiments are optimized for predicting the onset of generalized idiopathic 
and/or symptomatic epilepsies. Generalized epilepsies usually are not localized to a focal 

10 point in the patient's brain as the patient's whole brain is affected by ictal episodes. As a 
result, generalized idiopathic epilepsies often produce more generaUzed symptoms in the 
patient. Examples of generalized idiopathic and/or symptomatic epilepsies include, but are 
not limited to, benign myoclonic epilepsy in infants, juvenile myoclonic epilepsy, childhood 
absence epilepsy, and juvenile absence epilepsy, infantile spasms (West syndrome), 

15 Lennox-Gastaut syndrome, and progressive myoclonus epilepsies. 

In still further embodiments, the present invention provides systems and methods 
directed to predicting ictal onset in patients afflicted with unclassified and cryptogenic 
epilepsies such as febrile convulsions, Landau Kleffher syndrome, and Rasmussen's 
syndrome. 

20 In some embodiments, the present invention provides a system for predicting ictal 

onset in a subject comprising: a first data sensor positioned on the scalp of a subject near the 
focal point of ictal onset; a second data sensor positioned on the scalp of the subject, 
wherein the second data sensor is remote from the first data sensor; and a processor 
configured to analyze data collected from the first and the second data sensors to provide a 

25 nonlinear mathematical manipulation of the data collected from the first and from the 
second data sensors, wherein the nonlinear mathematical manipulation produces a first 
marginal predictabiUty value, and a second marginal predictability value. Preferably, the 
first and the second data sensors comprise electrodes. The present invention is not limited 
by the number, type, or placement of electrodes on a subject. In some embodiments, a 

30 plurality of electrodes (e.g., 2 ... 10 ... 100 . 200 ... or more) may be placed about (e.g., 
implanted or attached (e.g., glued, taped, etc.) to the subject's skin). It is further 
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contemplated that those skilled in the art of electrode placement can determine the suitable 
type, number, and placement of electrodes to best predict the seizures in a particular subject. 

In preferred embodiments, electrodes are used to collect electroencephalogram 
(EEG) data from subjects. The present invention is not limited however to collecting EEG 
5 data from subjects. Indeed, other types of data may be collected and used to help predict 
seizures in a subject. For example, some embodiments of the present invention contemplate 
using, MEG, and/or ECoG data. In still further embodiments, the methods and systems of 
the present invention are configured to incorporate subject data collected from other sources 
such as magnetic resonance imaging (MRI), x-rays, including, computed tomographic (CT) 

10 scans, various genotypic and phenotypic based tests and observations, including, but not 
limited to, comparative genomic hybridization, polymerase chain reaction based 
microsatellite analysis, fluorescence in situ hybridization studies, tissue biopsies, physical 
examinations, and the like. 

In preferred embodiments, the methods and systems of the present invention 

1 5 comprise a processor component (e.g., a computer comprising a computer readable memory 
and/or a devices to establish a communications link with one or more other component 
devices of the present systems) configured to compare/analyze the difference between a first 
marginal predictability value and a second, or more, marginal predictability value. In some 
embodiments, the marginal predictability values are computed by manipulation of subject 

20 data using nonlinear techniques and manipulations. However, the present invention is not 
limited to only using nonlinear data manipulation techniques. Indeed, in certain 
embodiments, one or more additional analysis techniques/values are considered, such as 
linear manipulations, calculation of various entropies, wavelets (i.e., portion(s) of the EEG 
with certain temporal and frequency characteristics) are analysis. 

25 In particularly preferred embodiments, difference determined between a first 

marginal predictability value and the second, or more, marginal predictability value 
decreases as a seizure, or a set of seizures, in a subject approach thus indicting ictal onset. 

Additionally, the present invention also provides methods and systems comprising 
one or more subject waming devices. In some of these embodiments, the subject warning 

30 devices are configured to receive data from a processor (e.g., computer processor) predicting 
that a seizure, or a set of seizures, in the subject is likely approaching. Accordingly, 
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preferred embodiments of the present invention provide warning to the subject (or a third 
party, such as, medical professionals, first responders, caretakers, and the like) that a 
seizure, or a set of seizures, is likely approaching. A number of subject warning devices are 
contemplated including, but not limited to, audible, visual, and tactile alarms. 
5 Some embodiments of the present invention further provide, systems comprising at 

least one component device in communication with a processor that administers anti-seizure 
agents to the subject. Anti-seizure agent administering devices suitable for use in the 
present invention include, but are not limited to, drug delivery devices (e.^., pumps, 
micropumps, patches, and the like), and electrical stimuli generators, and the like. In some 

10 embodiments, the component devices {e.g.y electrodes, component device for administering 
anti-seizure agents, processors, communications links, power supplies, etc.) are implanted in 
the subject; in some other embodiments, they are not implanted. 

In yet other embodiments, the methods and systems of the present invention provide 
additional component devices and commimications links to one or more third parties that 

15 transmit information from the system predictive of a seizure, or set of seizures, in the 
subject. 

In additional embodiments, the present invention provides methods for predicting 
ictal onset in a subject comprising providing: a subject; a system configured to detect ictal 
onset, wherein the system comprises a first data sensor positioned on the scalp of the subject 

20 near the focal point of ictal onset; a second data sensor positioned on the scalp of the 
subject, wherein the second data sensor is remote firom the first data sensor; a processor 
configured to analyze data collected firom the first and the second data sensors to provide a 
nonlinear mathematical manipulation of the data collected from the first and fi-om the 
second data sensors, wherein the nonlinear mathematical manipulation produces a first 

25 marginal predictability value, and a second marginal predictability value; and a subject 
warning device in communication with the processor; and contacting the subject with the 
system; determining the fist marginal predictability value and a second marginal 
predictability value; predicting ictal onset in the patient by difference in the fist marginal 
predictability value and a second marginal predictability value. 

30 
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DESCRIPTION OF THE FIGURES 

Fig. 1 shows a schematic representation of some of the standard scalp (Intemational 
10-20 System) and sphenoidal electrode coverage over the left hemisphere. Right 
hemisphere electrodes match those displayed, but have the next highest even number. 
5 Fig. 2A shows 62 as a function of time for a temporal electrode (F8) ipsilateral to the 

side of ictal onset during a 20 minute interictal period. 

Fig. 2B shows 82 as a function of time for the occipital (02) electrode ipsilateral to 
the side of ictal onset during the same period as in Fig. 2A. 

Fig. 3 A shows 62 as a function of time for a temporal electrode (F8) ipsilateral to the 
10 side of ictal onset during a one hour preictal period leading up to a seizure 

Fig. 3B shows 82 as a function of time for the occipital (02) electrode ipsilateral to 
the side of ictal onset during the same period as in Fig. 3 A. 

Fig. 4A shows Q2 as a function of time during the same period as in Fig. 2A. 

Fig. 4B shows Q2 as a function of time during the same period as in Fig. 3A. 
15 Fig. 5 A shows 82 as a function of time for a temporal electrode (F7) of a non- 

epileptic subject during a 20 minute period. 

Fig. 5B shows 8 2 as a function of time for the occipital electrode (Ol) of a non- 
epileptic subject during the same period as in Fig. SA. 

Fig. 5C shows Q2 as a function of time during the same period as in Fig. 5A. 
20 Fig. 6 shows the sum of positive ranks (SPR) for a collection of 40 interictal epochs 

of 20 minute duration taken from 8 patients with TLE. 

Fig. 7 shows the SPR test for 24 one hour preictal epochs prior to seizure onset from 
the same cohort of patients as in Fig. 6 

Fig. 8 shows the SPR for 24 epochs of 20 minute duration taken from 6 non- 
25 epileptic subjects. 

Fig. 9 shows Q2 for 4 interictal and 3 preictal epochs for one patient. 

DEFINITIONS 

To facilitate an understanding of the present invention, a number of terms and 
30 phrases are defined below. 

7 
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As used herein, the term "epileps/* refers to a heterogenous condition comprising 
recurrent unprovoked seizures characterized by abnormal transient paroxysmal disturbances 
ie.g.^ excitability) of cerebral cortical neurons, which leads to paroxysmal discharges in 
neuronal aggregates. Epileptic distiurbances in brain function can manifest as episodic 
5 impairment or loss of consciousness, abnormal motor phenomena, psychic or sensory 
disturbances or perturbation of the autonomic nervous system. On the basis of origin, 
epilepsy is characterized as being either idiopathic (Le., essential and/or genetic) or 
symptomatic acquired, secondary to a known insult such as a stroke or head trauma, or 
where the acquired cause is of uncertain origin cryptogenic]). On the basis of clinical 

10 and electroencephalographic observations, four common subdivisions of epilepsy are 

recognized. One, idiopathic generalized, epilepsies of unknown cause but usually with a 
heredofamilial basis. Examples include, but are not limited to, juvenile absence, pure 
absence, and juvenile myoclonic epilepsy. Two, symptomatic generalized, epilepsies due to 
a know insult such as birth hypoxia or of uncertain insult (commonly referred to as 

15 cryptogenic). Examples include, but are not limited to. West Syndrome, and Lennox 
Gastaut Syndrome. Three, idiopathic localization-related (synonymous with idiopathic 
partial or focal), epilepsies believed to be of local origin in a part of flie brain but likely due 
to heredofamilial causes. Examples include, but are not limited to, benign epilepsy of 
childhood with centrotemporal spikes (Benign Rolandic Epilepsy) and benign epilepsy of 

20 childhood with occipital spikes. Fovu", symptomatic localization-related (Symptomatic 

Partial or Focal), epilq)sies commonly believed to be of local origin in the brain and due to 
an identified or cryptogenic insult. Examples include, but are not limited to, temporal lobe 
epilepsy, frontal lobe epilepsy, occipital, or parietal lobe epilepsy. Temporal lobe epilepsy 
is commonly subdivided into those seizures arising from mesial temporal lobe structures, 

25 and those from neocortical parts of the temporal lobe. Other subdivisions of epilepsy 
arising in other areas is also possible. 

As used herein, the terms "seizure" or "ictal onset" refer to abnormal involuntary 
paroxysmal electrical activity (e.g., discharges of neurons) in the cerebral cortex that alters 
neurological function. Seizure can result in wide variety of clinical manifestations, 

30 including but not limited to, muscle twitches, staring, tongue biting, involuntary urination, 
loss of consciousness, and total body shaking. 
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As used herein, the term "subject" refers to an animal provided with the systems and 
methods of the present invention. Such organisms preferably include, but are not limited to, 
mammals (e,g,y murines, simians, equines, bovines, porcines, canines, felines, and the like) 
and most preferably humans. In the context of the invention, the term "subject" generally 
5 refers to an individual having a neurological condition (e.g., epilepsy) contemplated to be 
benefited by use of the systems and methods of the present invention for seizure prediction. 

The term "diagnosed," as used herein, refers to the recognition of a disease by its 
signs and symptoms (e.g.y uncontrollable neiu^ological activity), genetic analysis, 
pathological analysis, or histological analysis, and the like. 
10 As used herein, the term, "m v/Yro" refers to an artificial environment and to processes or 
reactions that occur within an artificial environment. In vitro environments can consist of, 
but are not limited to, test tubes and cell cultures. 

The term "«n vivo" refers to the natural environment (e.g., an animal or a cell) and to 
processes or reactions that occur within a natural environment. 
15 In preferred embodiments, the term **target cells" refers cells and tissues 

contemplated for use in the present invention. "Target cells" include, but are not limited to, 
neural cells. 

As used herein, the term "marginal predictability** refers to the quantity defined in 
formula 6, below. It is a measure of the extent to which the dth lag in a time series or other 

20 data set provides additional predictive information for the next term in the series, give that 
the information in the intervening d-1 lags has already been used. 

As used herein, the term "electroencephalograph (EEG)," and grammatical 
equivalents, relates to systems and devices for recording electric potentials of the brain 
derived from electrodes attached to a subject's scalp. Electroencephalographic equipment 

25 and techniques are used for collecting electroencephalographic data from a subject (e.g., a 
patient afflicted with epilepsy). This data is used by the systems and methods of the present 
invention to predict the onset of seizures (e.g., ictal onset) in the subject. The American 
Clinical Neurophysiology Society (ACNS, Bloomfield, CT) provides technical guidelines 
for performing EEGs. 
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The term "data sensors," As used herein, refers to scalp based electrodes used for 
collecting data {e,g, , EEG, EMG, and the like) from the subject that is subsequently used by 
the systems and methods of the present invention for predicting ictal onset in the subject. 

As used herein, the term "intracranial data sensors," refers to electrodes that are 
5 implanted in a subject's cranium used for collecting data (e.g., EEG, EMG, and the like) 
from the subject that is subsequently used by the systems and methods of the present 
invention for predicting ictal onset in the subject. 

As used herein, the term "subdural data sensors," refers to electrodes that are 
implanted under a subject's dura mater used for collecting data (e.g., EEG, EMG, and the 
10 like) from the subject that is subsequently used by the systems and methods of the present 
invention for predicting ictal onset in the subject. 

As used herein, the terms "subject warning device," or "patient warning device" 
refer to devices configured to receive information from a processor indicative of ictal onset 
and then to use this information to activate features of that device that alert a subject to the 
1 5 likelihood of ictal onset (e.g,, a seizure). "Subject warning devices" suitable for use in the 
present invention include, but are not limited to, audible (e,g,y buzzers, beepers, and the 
like), tactile (eg., vibrating, pulsating, and the like), and visual alarms (e.g., flashing and the 
like). 

As used herein, the term **preictal period" refers to the period of time preceding the 
20 onset of seizures (e.g., ictal onset) in a subject calculated by analysis of subject data (e.g., 

EEG, medical examination) that indicates that the subject is likely to start seizing in the near 

ftiture (e.g., 1-5, 10, 30, 60, 120 min. or more). Likewise, the term "preictal warning 

period" refers to a portion (or all) of the "preictal period" where the subject is alerted to the 

likelihood of an ictal onset. 
25 As used herein, the term "interictal" refers to the period of time between the 

termination of an isolated seizure or a set of seizures and the beginning of a second isolated 

seizure or set of seizures (e.g., ictal onset). 

As used herein, the term "ictal onset" refers to the time or short epoch at which a 

seiziu-e begins in a subject. The period of "ictal onset" is usually preceded by a "preictal 
30 period" of sufficient duration such that the systems and methods of the present invention 

provide the subject with a sufficient "preictal warning period" to begin administration of 

10 
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anti-seizure therapies and/or to prepare (e.g., stop/abstain from operating machinery etc.) for 
the onset of seizures. 

As used herein, the temi "focal point of ictal onset" refers to location(s) in a 
subject's brain where abnormal discharges (e.g., seizures) originate and spread forth from. 
5 As used herein, the term "post ictal" refers to the period of time immediately 

following the cessation of seizures in a subject to a point in time when the subject has 
returned to baseline levels of neurological frmction. 

As used herein, the term, "anti-seizure agent" refers to a compound or to an 
electrical stimulus used to treat epilepsy in a subject. Treatment is meant to encompass 
10 preventing, delaying, or lessening the severity of seizures and generally controlling epilepsy 
and its effects in a subject. 

As used herein, the term "anti-seizure agent administering device," refers to a device 
configured to administer at least one anti-seizure agent to a subject upon receiving 
instruction (e.g., via sl communications link and/or user interface, etc.) provided by the 
15 systems of the present invention (e.g., a computer processor), the subject, a healthcare 
provider, a first responder, or a caretaker, etc. 

As used herein, the term "effective amount" refers to the amoimt of an agent (e.g., an 
anti-seizure dmg) or to the characteristics of the voltage of electrical stimuli sufficient to 
produce beneficial or desired effects in a subject (e.g., a subject having epilepsy). An 
20 effective amount of a dmg or other therapeutic can be provided in one or more 

administrations, applications, or doses and is not intended to be limited to particular 
formulation or administration routes. 

As used herein, the term "modulate" refers to the activity of an anti-seizure agent 
such as a medicament (e.g., an anti-seizure drug) or other therapeutic agent (e.g., electrical 
25 stimuli) to affect (e.g. , to promote or retard) any aspect of neural cell function, including, 
but not limited to, electrochemical discharge and the like. 

As used herein, the term "program," when used as a noun, refers to an organized list 
of instructions that, when executed performed by the computer) causes the computer to 
behave in a predetermined manner (e.g., to calculate and compare the difference in marginal 
30 predictability values for each of the two time series). Programs consist of modules, each of 
which contains one or more "routines/subroutines" (i.e., a section of a program that 

11 
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performs a particular task). In preferred embodiments, the programs described herein are 
written in a high-level programming language, such as BASIC, FORTRAN, C, C++, 
PASCAL, COBOL, or LISP, and the like. 

As used herein, the terms "processor," "computer processor," "computer," and 
5 "central processing unit," or "CPU" are used interchangeably to refer to a device that is able 
to read a program from a computer memory (e.^., ROM or other computer memory) and 
perform a set of steps according to the program. The systems and methods of the present 
invention for seizure prediction are not intended to be limited, however, to implementation 
by processors. Indeed, in some embodiments the present invention contemplates systems 

10 comprising logic chips, logic arrays, gate arrays, application-specific integrated circuits, and 
programmable logic devices (PLDs), and the like. 

As used herein, the terms "memory device", "computer memory device", or 
"computer memory" refer to any storage media readable by a computer processor. 
Examples of computer memory include, but are not limited to, RAM, ROM, computer 

15 chips, digital video disc (DVDs), compact discs (CDs), hard disk drives (HDDs), and 
magnetic tape. 

As used herein, the term "computer readable medium" refers to any device or system 
for storing and providing information (e.g., data and instmctions) to a computer processor. 
Examples of computer readable media include, but are not limited to, DVDs, CDs, HDDs, 

20 and magnetic tapes. 

The term "operably linked," in one sense, when used in reference to the operation of 
the disclosed systems and methods for predicting ictal onset refers to the execution (e.g., 
performance by a computer) of a computer program by a computer processor (or other 
suitable devices such as, logic chips and the like) to produce a desired result (e.g., the 

25 prediction of ictal onset in a subject). In another sense, the term "operably linked," when 
used in reference to the operation of the disclosed systems and methods, refers to computer 
hardware devices and other apparatuses such as, communications links [e.g., fiber optics, 
modems, infrared [IR], LANs, the Internet]) configured to receive and/or exchange 
information with the disclosed systems, methods and programs for predicting ictal onset. In 

30 preferred embodiments, operably linked computer hardware devices and other apparatuses 
are configured to receive and/or exchange information with a computer program stored in 

12 
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computer readable memory associated with a computer processor via, for example, wires or 
cables, computer cards and boards, circuits, communication links fiber optics, 
modems, IR, LANs, the Internet], etc.^ and any necessary device drives or computer 
program subroutines stored in computer readable memory. 

5 As used herein, the terms "user," and "system user" when used in reference to 

controlling the operation of a computer program (e.g., a software program stored in 
computer memory that predicts ictal onset), refers to a person, or a second computer 
program and system (e.g., a software program stored in computer memory), that controls the 
operation of the first computer program by selecting and/or entering system operation 

10 parameters and information. In some embodiments, the **user" is also the "subject." 

As used herein, the term "user interface" refers to the junction between a user and a 
computer program (e.g., the user interface is configured to be capable of receiving 
information fi-om the user). A user interface allows the user to transmit or convey 
commands (e.g., instructions) to a computer program, hardware device or apparatus to 

15 perform specific tasks (e.g., prediction of ictal onset), "Graphical user interface," or (GUI), 
refers to a user interface that takes advantage of the computer's graphics capabilities for 
entering and retrieving data from a program. 

The term "Internet", as used herein, refers to a collection of interconnected (public 
and/or private) networks that are linked together by a set of standard protocols (such as 

20 TCP/IP and HTTP) to form a global, distributed network. While this term is intended to 
refer to what is now commonly known as the Internet, it is also intended to encompass 
variations that may be made in the future, including changes and additions to existing 
standard protocols. 

As used herein, the terms "World Wide Web" or "Web" refer generally to both (i) a 
25 distributed collection of interlinked, user-viewable hypertext documents (commonly referred 
to as Web documents or Web pages) that are accessible via the Internet, and (ii) the client 
and server software components which provide user access to such documents using 
standardized Internet protocols. Currently, the primary standard protocol for allowing 
applications to locate and acquire Web documents is HTTP, and the Web pages are encoded 
30 using HTML. However, the terms "Web" and "World Wide Web" are intended to 
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encompass future markup languages and transport protocols that may be used in place of (or 
in addition to) HTML and HTTP. 

DESCRIPTION OF THE INVENTION 

5 Epilepsy is often though of as being a symptom of some xmderlying brain 

dysfunction that causes the dysfunction of brain cells, either at a specific place (as in focal 
epilepsy), or more widely (as in generalized epilepsy). This neurological dysfunction causes 
physical consequences such as the loss of consciousness or loss of muscle control. 

Despite the difficulties in classifying different types of epilepsies the medical 

10 sciences have made several useful distinctions in the types of epilepsies. For example, 

partial epilepsies are characterized as originating from a clearly defined focal area within the 
brain. As a result, partial epilepsies have symptoms characteristic of their site of origin, 
such as simple visual hallucinations (for seizures of the occipital lobe) and unilateral motor 
difficulties (for seizures arising from the frontal lobe). 

15 The present invention relates to methods and systems/devices for noninvasive 

nonlinear prediction of ictal onset in patients afflicted by neurological conditions. In 
particular, the present invention provides methods and systems for noninvasive nonlinear 
prediction of seizures in patients afflicted with epilepsy. In particularly preferred 
embodiments, the systems and methods of the present invention are directed to predicting 

20 seizures (e.g., ictal onset) in a subject afflicted with a focal (localization-related or partial) 
epilepsy. 

The systems and methods preferably being based on analysis of two or more 
electroencephalogram (EEG) recordings, one set of recordings taken from an electrode close 
to the region of ictal onset, and a second or more set of recordings (e.g.y concurrent 

25 recordings) taken from a region remote from the region of ictal onset. The respective 
recordmgs of subject data (e.g., EEG data) are each converted into time series. Either or 
both of the electrode locations may be ipsilateral to (i.e., on the same side as), or 
contralateral to (/.e., on the opposite side as) the focal point of seizure onset. In particularly 
preferred embodiments of the present invention, a comparison is made between the EEG 

30 data retrieved by a first electrode adjacent to the site of seizure onset (e.g., temporal 

location) and second electrode remote (e.g.y occipital location) from the site of seizure onset. 
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However, in some other embodiments of the present invention, different combinations and 
placements of electrodes are contemplated. In some embodiments, electrodes are placed 
about a subject following the standard 10-20 placement system. However, those skilled in 
the art appreciate that alternative electrode placements are preferable in recording data from 
5 certain types of epilepsy or other seizure causing disorders. Indeed, while collecting and 
manipulation of a subject^s EEG data is preferable in some embodiments, collection and 
manipulation of other types of data are contemplated in additional embodiments {e.g,, MEG, 
ECoG, and the like). 

A nonlinear quantity called the marginal predictability (MP) is calculated for each of 

10 the time series generated using the formulas described herein. These MP values are then 
compared. It is contemplated that the difference between the two MP values predictably 
decreases several tens of minutes (e.g,, 10, 20, 60 min., or more) prior to ictal onset. It is 
during a portion of the preictal period that the subject is preferably wamed of the potential 
for an oncoming seizure (ictal onset) by one or more subject warning devices. In still other 

15 embodiments, the MP values (or other related quantities) are calculated based on 
intracranial or subdural data recordings. 

The systems and methods of the present invention are not intended to be limited to 
incorporating any particular type of subject warning device(s). For example, in some 
embodiments, suitable warning devices for use with the seizure prediction systems and 

20 methods include, but are not limited to, audible (e.g., beeping devices), tactile (eg., 

vibrating devices), visual warning devices (e.g., flashing devices). The present invention 
contemplates that one or more of these types (or other types) of warning device be 
incorporated into the systems and methods of the invention. In some embodiments, the 
subject waming devices provide the subject with an escalating waming signal (e.g., 

25 increasing volume audible alarm) as the onset of an ictal episode draws closer (e.g., 
temporally near). 

In preferred embodiments, subject data is collected on one or more channels (e.g., 1, 
2 ... 4 ... 8 ... 16 ... 32 .. . 128 . . . 256 ... or more). In preferred embodiments, the foci 
of a subject's seizures are mapped, for example, using the International 10-20 system of 
30 electrode placement. In preferred embodiments, the number and placement electrodes (data 
sensors) used with a particular subject correspond to the determination of seizure foci. In 
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particularly preferred embodiments, the systems and methods of the present invention 
collect patient data from two channels. 

Some embodiments of the present invention use electrodes specially designed to 
reduce noise acquired from the data acquisition process (e.g., system noise, subject muscle 
5 contractions, etc). For example, in one embodiment, the systems use one or more 

electrodes that integrate a first amplifier stage with a sintered Ag-AgCl electrode to provide 
extremely low-noise and virtually interference free extracranial data sensing (e.g., BioSemi, 
Amsterdam, Netherlands). In preferred embodiments, extracranial electrodes are used with 
noise levels near the thermal noise level of electrode impedance. 

1 0 Preferably, the other components of the systems of the present invention are selected 

based on considerations, including, but not limited to, size (e.g., miniaturization), weight, 
power consumption, impedance, signal-to-noise ratio and the like. 

In still fiirther embodiments, subject data is transmitted via a communications link 
(e.g., telephone line, wireless network (e.g., radio wave, infra red, microwave, etc), coaxial 

15 cable, fiber optic cable, wire, circuit path, and the like) to a processor configured to analyze 
the subject data and to determine the temporal proximity of an ictal episode. The 
communication links of the present invention are not intended, however, to be Umited to any 
particular frequencies or wavelengths. Thus, in some embodiments, the processor is located 
remote from the subject. In other embodiments, subject data acquisition device(s) (e.g., 

20 EEG electrodes) are integral to the processor used. In preferred embodiments, all of the 
devices/systems and methods of the present invention are integrated and sufficiently 
miniaturized such that they are easily carried by the subject. 

In additional preferred embodiments, the systems and methods of the present 
invention further incorporate data storage devices. It is contemplated that the data storage 

25 devices be configured to store the programs for predicting seizures from subject data. It is 
fiirther contemplated that the data storage devices be configured to record subject data (e.g., 
EEG recordings) such that the data is transferable (e.g., via a communications link) to one 
or more additional devices. Data storage devices suitable for use in the systems of the 
present invention include, but are not limited to, floppy disks, hard disks, random access 

30 memory, digital tapes, compact disks, digital video disks, and the like. In preferred 

embodiments, subject data is periodically (e.g., in real time, hourly, daily, weekly, monthly. 
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etc) transferred to the subject's healthcare provider or to another third party (e.g., 
emergency responders, insurers, employers, and the like). As described above, subject data 
is transferred by any suitable communications link including encrypted or unencrypted 
transmission over the Internet. Thus, in certain embodiments, the systems and programs of 
5 the present invention for predicating ictal onset are accessible via the Internet (z.e., World 
Wide Web) over a communication network. Prior to transmission, transfer, or storage 
subject data is preferably compressed using one or more software compression programs. 

Certain embodiments of the systems of the present invention further provide one or 
more component devices for administering therapeutic anti-seizure (e.g., dmgs and/or 

10 electric stimuli) agents/treatments to the subject. In preferred embodiments, administrations 
of anti-seizure agents/treatments are timed to coincide with the system's prediction of ictal 
onset. In other embodiments, administrations of anti-seizure agents/treatments are at regular 
scheduled intervals (e.g., hourly, daily, weekly, etc.). Of course, additional variations on the 
administration times and routes of anti-seizure agents/treatments are within the scope of the 

15 present invention. For example, in some additional embodiments, component devices are 
configured to administer anti-seizure agents on demand. In still other embodiments, a 
conmiand to the anti-seizure agent administering device instructing it to administer the anti- 
seizure agent is sent from a third device (e.^., a processor) or entity remote from the seizure 
prediction systems via a commimications link (e.g., telemetry link). In yet other 

20 embodiments, the seizure prediction systems of the present invention optionally incorporate 
additional programs and devices for determining the proper time(s) to administer anti- 
seizure agents to increase their effectiveness. 

The present invention further contemplates a number of component devices that re 
suitable for administering anti-seizure agents to subjects. For example, in some 

25 embodiments, the systems of the present invention incorporate electric stimulus devices 
(e.g., microelectromechanical systems (MEMS), or wireless integrated microsystems 
devices (WIMS), and the like) configured to provide the subject with mild electrical pulses 
prior to ictal onset. For example, in some embodiments, the present systems incorporate 
devices configured to provide the subject with adaptive electric field therapy or vagus nerve 

30 stimulation (VNS) therapy (Cyberonics, Inc., Houston, TX). In some embodiments, the 
component devices for delivering electrical stimuli are implanted in the subject. In other 
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embodiments, component devices for delivering electrical stimuli are not implanted in the 
subject. In particularly preferred embodiments, these devices are connected via a 
communications link the seizure prediction systems of the present invention. Optionally, 
these devices are also connected via a communications link to remote processors controlled 
5 by the subject's healthcare provider and the like. 

Those skilled in the art appreciate that a number of implantable medical devices 
capable of delivering electric stimuli to a subject are suitable for use in the present systems 
(e.^., micro-electro-mechanical systems (MEMS)). The present invention is specifically not 
intended, however, to be limited to incorporation of implantable devices. 

10 A number of systems for powering and controlling implantable/non-implantable 

medical devices are suitable for in the present invention. For example, some implantable 
component devices suitable for use in the present invention incorporate wireless power 
transfer systems that operate through an inductive link to a battery {e,g, , rechzirgeable 
Hthium battery) using bi-directional data communication. In other embodiments, 

15 implantable components such as micro pumps (actuators) and/or electric stimulators are 
configured to run on integral batteries. In particularly preferred embodiments, the batteries 
are energized using an RF antenna or low-fi-equency magnetic loop implanted in the 
subcutaneous tissue. 

In still fiirther preferred embodiments, the components of the present invention 

20 implanted in a subject are hermetically sealed and constructed of non-immunogenic 
materials. 

Anti-seizure drugs suitable for use in some embodiments of the present invention 
include, but are not limited to, phenobarbitol, phenytoin, carbamazepine, ethosuximide, 
valproate, benzodiazepines (e.g., clonazepam, diazepam, and lorazepam), felbamate, 

25 gabapentin, lamotrigine, vigabatrin, topiramate, and the like. In preferred embodiments, one 
or more anti-seizure drugs are administered to a patient by a component device under the 
instruction of the systems of the present invention. For example, in a preferred 
embodiment, when the systems of the present invention detect ictal onset is likely, the 
system instmcts (e.g., via a communications Hnk) a component device (e.g., a micropump) 

30 to administer a therapeutically effective amount of an anti-seizure drug to the subject. In 
other embodiments, the systems of the present invention periodically (or aperiodically) 
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instructs a component device (e.g., micropump) to administer a therapeutically effective 
amount of an anti-seizure drug based on a preset schedule (e.g,, hourly, daily, weekly, etc) 
or upon the subject's self instruction, or upon being instructed to do so by a third party (e.g., 
the subjects health care provider), or upon being instructed to do so by the invention which 
5 may instruct intervention in an aperiodic, or therapeutically random fashion. In some of 
these embodiments, the third party health care provider is remote from the subject and 
instructs a component device to administer an anti-seizure dmg to the subject via a 
communications link. 

In particularly preferred embodiments, the present invention provides methods and 

10 systems for predicting the onset of seizures in a subject sufTering from a focal epilepsy (e.g., 
mesio-basal temporal lobe epilepsy). However, additional embodiments are provided for 
predicting the onset of seizures in subjects with various types of non- focal epilepsies. 

In preferred embodiments, subject data tapes are collected from scalp-sphenoidal 
monitoring (e.g., data sensors) of epilepsy patients. Preferably, the subject data (e.g., EEG) 

15 is sampled at 200 Hz and the data tapes are read on a Viewmaster 5000 machine 

(BMSI/Nicolet, Madison, WI). The subject data is transferred via a communications link to 
a computer memory for storage. From the computer memory the data is transferred by via a 
communications link (e.g., ftp) to a computer workstation (e.g., Linux workstation). The 
subject data is compressed, such that each data point takes up about 2 bytes. In some 

20 embodiments, for purposes of analysis, the compressed subject data is extracted and 

converted from binary to ASCII using standard processes. Further data analysis is carried 
out on the computer workstation. In some embodiments, the subject data set used in the 
analyses of the present invention are decimated, and consist of every third data point from 
the original 200 Hz sampled data record. However, the decimation frequency is a parameter 

25 that can be changed in the analysis. Additional filtering of the original subject data is 
typically not necessary, however, in certain other embodiments original subject data is 
subjected to additional or different filtering steps. 

Preferred embodiments of the systems and methods of the present invention 
comprise ambulatory devices in which one or more data sensors (e.g., electrode leads) are 

30 attached to the subject's scalp. The data collected from the several leads (e.g., EEG data) is 
temporarily stored for form one to several, or more, minutes in one or more temporary 
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Storage devices {e.g,, computer memory) worn by the subject. These devices also contain a 
processor to compute the nonlinear and other quantities of interest. In some embodiments, 
once the data has been used for the calculations of interest, that data is purged and replaced 
.with new subject data coming in real time from the subject's data sensors (e.g., electrodes). 
5 In other embodiments, the data is stored for longer periods of time and used in subsequent 
analyses in one or more temporary devices worn by the subject. The present invention is not 
limited however to storage devices, or other components, of the present invention that are 
wom by the subject. In still other embodiments, the subject data is retained and 
subsequently transferred via a communications link to additional storage devices and 

10 processors. In some of these embodiments, the subject data is optionally further compressed 
for easier transfer and storage. The computed value of the nonlinear (and possibly linear) 
quantities of interest are added to an existing record and are used with recently 
analyzed/manipulated (e.g., about the past 1 ... 5 ... 10 ... 30 ... or 60 min. or more) 
subject data to generate seizure predictions. In preferred embodiments, the time necessary 

15 to perform these calculations is such that new values of the nonlinear metrics are computed 
with only a few (e.g., 1 ... 5 ... 10 ... to about 20) seconds lag from real time. 

In some embodiments, the seizure prediction systems are attached physically, 
electronically, or by wireless communication to other devices(s) that provide medical or 
electrical intervention to abort a seizure or to ameliorate its effects. 

20 Additional exemplary systems and methods of the present invention are described in 

more detail in the following sections: 1. Mathematical considerations; and n. Data analysis. 

L Mathematical considerations 

The first attempts at seizure prediction relied on standard linear statistical methods. 

25 (See e.g., J. Gotman et al.. Long-term monitoring in epilepsy, EEG suppl. No. 37 

Amsterdam: Elsevier [1985], and 1. Osorio et al. Epilepsia, 39:615-627 [1998]). Recent 
efforts by brain researchers have included finding nonlinearity and claims of chaos in brain 
study data sets (e.g., EEGs). (See J. Theiler, Phys. Lett. A., 196:335-341 [1995]; A. 
Babloyantz and A. Destexhe, Proc. Natl. Acad. Sci. U.S,A., 83:3513-3517 [1986]; G.W. 

30 Frank et al, Physica D, 46:427-438 [1990]). However, many attempts at applying nonlinear 
dynamics and chaos theory to biological systems have met with limited success. (See e.g, 
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Theiler supra; J. Theiler and P. Rapp, Electroenceph. Clin. Neurophysiol., 98:213-222 
[1996]; M. Palus, Nonlinearity in normal human EEG: Cycles, temporal asymmetry, non- 
stationarity and randomness, not chaos, Santa Fe Institute working paper 94-10-054 [1994]); 
Babloyantz and Destexhe supra, and L. D. lasemidis et al. Brain Topog. 2:187 [1990]; A. 
5 Destexhe and A. Babloyantz, Deterministic chaos in a model of the thalamo-cortical 
system. In: Self-Organization, Emerging Properties and Learning, A. Babloyantz ed.. 
Plenum Press, New York, NY [1991]). Thus, prior to the present invention effective 
methods and systems for predicting ictal onset in subjects using noninvasive nonlinear 
analyses were imknown. For instance, most of the research concerning prediction of 

10 epileptic seizures has focused on collecting and analyzing intracranial data recordings. 
Intracranial data recordings have traditionally been available from a subset of epilepsy 
patients with medically refractory temporal lobe epilepsy (TLE). In contrast, however, 
preferred embodiments of the present invention collect and analyze subject data obtained 
from extracranial (e.g., scalp) electrodes. 

15 However, in some other embodiments, the present invention contemplates placement 

of intracranial electrodes (e.g., in cases where the seizure focus is known). 

In preferred embodiments of the present invention, nonlinear mathematical methods 
are used to manipulate subject data to predict the onset of seizures in an epilipetic subject. 
The basis for the nonlinear dynamical analyses used in the present invention is the 

20 correlation integral (P. Grassberger and 1. Procaccia, Physica D, 9:189-208 [1983]) defined 
as: 

CAyii)^yUy) = P(\y'''ii)-y''\j)\<^) (Formula 1) 

25 where P(0 denotes the probability of the argument, xj is the j^*" element of the time series 

being reconstructed, and y^'^^i) = (xi, Xi.|, . , ., Xi.d+i) is a d-dimensional vector reconstructed 
from data. The notation || " || means norm of the argument. In preferred embodiments, the 
methods of the present invention use the max norm equation, which is the computationally 
simplest definition, i.e., \y^''\i)-y^'^\j)\<z ifmax[|xi.k-Xj.k|] <8fork=0, 1, ...,d-l. The 

30 present invention is not intended, however, to be limited to using this defmition of norm. 
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Indeed, a number of different definitions can be used for the norm. The skilled artisan will 
recognize that good theoretical and/or computational reasons for choosing one definition 
over the other exist in particular embodiments. Thus additional embodiments of the present 
invention use other definitions for norm. The quantity Cd is the probability that two vectors 
5 reconstructed from the time series in d-dimensions will be close to each other. In terms of 
the original time series, Cd is a measure of the likelihood that tsvo sequences of length d 
taken from a time series will look similar. Using the Ca's, predictability can be defined as 
{See, R. Savit and M. Green, Physica D, 50:95 [1991]): 

10 5^=.^iL. (Formula 2) 

In view of formula 1 , Sd is the conditional probability equation 

= ^(^rf+i I 2, ) (Formula 3)' 

where, 

=|Xi+k.i - |<8 . (Formula 4) 

15 Thus, in preferred embodiments, Sd is the conditional probability that if two randomly 
chosen d-tuples from the time series have their first d-1 elements within 5 of each other, 
respectively, then the d* elements will also be within 5. 

In some embodiments, Sd is used as a nonlinear statistic. However, in preferred 
embodiments, a more sensitive discriminator of nonlinear structure in time series {See, R. 
20 Manuca et al. Mathematical Biosciences 147:1 [1996]); and R. Savit and M. Green, 
Physica D, 50:95 [1991]) is the ratio of Sd's, defined as: 



= = . (Formula 5) 



To make the interpretation simple, preferred embodiments of the present invention define 
25 marginal predictability as: 



Rd -1 

SrfS-^ — (Formulae) 

^d 
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Wherein, 8d is a measure of how much additional predictive information there is in the 
(d+1)*^ lag of the time series, given that information in the intervening d lags has already 
been used. It is contemplated that if 8d is close to zero, there is no additional predictive 
information, on average, for the current value of the time series in the value of the (d+l/^ 

5 lag. However, if 5d is significantly different fi-om zero, then Sd>Sd.i, and there is additional 
predictive information in the (d+1)^^ lag. "Predictive information** here is understood in the 
sense of nonlinear dynamics. (See e.g., R. Savit and M. Green supra; and K, Wu et al. , 
PhysicaD, 69:172 [1993]). 

In particularly preferred embodiments, the methods of the present invention 

10 comprise comparing 5d for two different scalp electrodes as a function of time. Thus, 

consider Qd(A,B;t) = 6d (A;t) - 5d(B;t), where A and B are two electrodes and t is time. In 
preferred embodiments, A will be one or more electrodes {e.g,, scalp electrodes) near the 
seizure focus and B will be an electrode {e.g., scalp electrodes) remote fi-om and preferably 
ipsilateral to the site of ictal onset. In the case of TLE, B will generally be an occipital 

15 electrode. Thus, preferred embodiments of the present invention compare differences in the 
marginal predictabilities of temporal and occipital electrodes between times far removed 
fix>m a seizure and times close to a seizure. In furdier embodiments, of the present 
invention the methods are based on similar calculations where i) A and B are contralateral 
temporal and occipital electrodes, respectively, and ii) A and B are temporal ipsilateral and 

20 contralateral electrodes, respectively. 

In analyses based on non-linear dynamics, several parameters {e.g, embedding 
dimension, and tolerance) need to be set prior to performing the calculations for determining 
the MP values. Thus, in some embodiments, the basic seizure prediction methods of the 
present invention can be tailored to the needs of a particular patient or for predicting 

25 particular types of epileptic seizures by adjusting the values of the parameters {e.g.^ altering 
the dimension of the reconstruction space) of the nonlinear mathematical manipulations. 

The present invention further contemplates additional embodiments wherein the 
methods of seizure prediction use one or more alternative mathematical quantities to 
increase the predictive value and robustness of the basic seizure prediction method in certain 

30 patients. For example, in some embodiments calculation of various entropies are 

contemplated. In other embodiments, wavelets (i.e., portion(s) of the EEG with certain 
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temporal and frequency characteristics) are analyzed to further predict seizures. In still 
other embodiments, additional linear or nonlinear quantities may be used alone or in 
conjunction with MP to increase the selectivity and sensitivity of seizure prediction. 

5 IL Data analysis 

In preferred embodiments, patient results are presented in 40 second windows 
comprising about 8000 data points (sampling rate of 200 Hz) and are used to calculate 6d 
and Qci(A.B;t). Preferably, the time series Xi is obtained from the originally sampled EEG 
recordings by using every third data point (t=3). This decimation of the data set was chosen 

10 to minimize the value of the mutual information. (See, H. Tong, Nonlinear Time Series 

Analysis: A Dynamical Systems Approach, Oxford University Press [1990]). e in (Formula 
4) is 10% of the standard deviation of the time series, x\. In other embodiments, it is 
contemplated that these values of x and e can be altered to improve predictability 
performance. Those skilled in the art appreciate the steps to alter this value without more 

15 than routine experimentation. In the results presented below, d=2. The results include 

twenty-four 20 minute preictal epochs from eight patients and forty interictal epochs of 20 
minute duration from eight patients in the analysis. Behavior states of the subjects were 
noted and epochs were chosen that represented a range of behavior states including various 
stages of awake alertness and various sleep states. Results of the analysis were independent 

20 of the behavior state of the subject. For purposes of comparison, twenty four epochs of 20 
minute duration from six non-epileptic subjects using methods of the present invention were 
also analyzed. 

Fig. 1 shows a schematic diagram indicating the standard nomenclature and 
placement of scalp electrodes according to the International 10-20 System. In particular, 

25 Fig. 1 shows the left side of a hypothetical patient's head. Electrodes placed in homologous 
locations on the right side of the scalp are labeled with the next highest even number. Thus, 
for example, the electrode on the right side of the head in the position homologous to the F7 
electrode is labeled F8. Depending on the side of the seizure focus, the F7 or F8 electrode is 
typically close to the site of ictal onset in cases of TLB. Two examples of 82 as a fimction of 

30 time are shown in Figs. 2A-B and 3A-B, respectively. Fig. 2A shows the time course of 82 
for a temporal electrode (F8) ipsilateral to the side of ictal onset during a 20 minute 
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interictal period. Fig. 2B shows 82 for the occipital (02) electrode ipsilateral to the side of 
ictal onset dviring the same period, 62(F8) is significantly higher than the 52(02). Fig. 3 A 
shows 82 for the same electrode (F8) used in Fig. 2A, but calculated during a one hour 
preictal period leading up to a seizure, and Fig. 3B shows 82 for the electrode (02) used in 
5 Fig. 2B, calculated during the same one hour preictal period used for Fig. 2A. Note that 
82(F8) is still greater than 82(02) until about 15 minutes before the seizure onset. Within 15 
minutes prior to the seizure, however, 82(F8) decreases to approximately the same level as 
82(02). 

Certain embodiments of the present invention are illustrated by considering Q2, 
10 which is the difference between the 82*8 of different electrodes. The results are presented in 
Figs. 4A-B. Note that in Fig. 4A, Q2 is greater than zero for the interictal epoch (typically 
close to 0.05). For most of the early portion of the preictal epoch Q2 is also greater than 
zero in Fig. 4B, but moves close to and stays near zero starting about 15 minutes prior to the 
seizure. Although there are differences in the profiles of the 82 's for different epochs, the 
15 features illustrated in Figs. 4A an 4B, namely, the fact that Q2 is smaller in the preictal 

compared to the interictal period can be found in almost all of the epochs studied from the 
epileptic subjects. 

The same analysis was also applied to a set of epochs taken from non-epileptic 
subjects. Figs. 5 A and 5C are examples of a 20 minute epoch from a non-epileptic subject. 
20 hi some of these embodiments, the electrodes comprise a temporal electrode (F7) and an 

occipital electrode (01), both on the left hemisphere. Fig. 5C shows that the value of Q2 for 
non-epileptic subjects are typically close to zero, which suggests that there is no systematic 
difference between the 82 of temporal electrodes and that of occipital electrodes for non- 
epileptic subjects. 

25 In order to statistically validate these observations, a Wilcoxon's sum of signed rank 

test was applied to the values of Q2. However, the present invention is not Hmited to the 
Wilcoxon's sum of signed rank test. Those skilled in the art recognize that other 
nonparametric tests for manipulating paired samples are also useful in certain embodiments 
of the present invention. Briefly, the Wilcoxon test (See e.g., E.L. Lehmann, 

30 Nonparametrics: Statistical Methods Based on Ranks, Holden-Day, Inc., Oakland, CA 
[1975]) is a nonparametric test for paired samples (Xi, Yi). Accordingly, in the present 
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methods the paired samples are the marginal predictabilities for focal electrodes and remote 
electrodes, respectively, i.e., (Xj, Yj) = (82(F electrodes), 82(0 electrodes). The Wilcoxon 
sum of signed rank test can be used to test the null hypothesis that the median of the 
difference, Dj = Xj - Yi, is equal to zero, so that it is just as likely that Xi>Yi as that Xi<Yi. 
5 Specifically, the test is calculated in the following way: 1) rank order the absolute values of 
the Di from smallest to largest (wherein Ri is the rank of | Di |, for i = 1, . . ., n); 2) assign the 
sign of Di to the rank of Dj; and 3) calculate SPR, (the sum of all the positive ranks), that is 
the sum of all those ranks that are associated with a positive value of Dj. If the above null 
hypothesis, namely that the median of the differences, D,, is zero, is true, about half the Dj 

10 values are positive and half negative, and SPR will be neither tool large nor too small, being 
close to n(n+l)/4, where n is the sample size. A test statistic can therefore be developed 
based upon SPR. Under the null hypothesis, the expected value of SPR is equal to n(n+l)/4, 
which is indicated by the broken lines in the Figs. 6, 7, and 8. In preferred embodiments, 
the probability of each distinct value of SPR under the null hypothesis may also be 

15 calculated, giving significance levels. For example, for n = 24, the probability that SPR is 
less than 81 under the null hypothesis is approximately 0.025. Hence, if SPR is less than 81, 
then rejection of the null hypothesis can be made with 97.5% confidence. 

The advantage of nonparametric test (e.g., the Wilcoxon signed rank test) over 
parametric tests (e.g., t-test) is that it does not make any ancillary assumptions about the 

20 distribution of Di. The only necessary assumption is that all Di's are indq)endently sampled 
firom the same distribution. 

Specifically, summed-rank tests were performed to test the null hypothesis that the 
median of 82 of the electrode adjacent to seizure onset, jiadi (e.g., electrodes F7 or F8, 
depending on which side of the brain contained the site of ictal onset), is the same as that of 

25 the electrode remote to the site of seizure onset Mremote (e.g., occipital electrode Ol or 02 
ipsilateral to the site of ictal onset), thus Q2=0, statistically. The simi of positive ranks 
(SPR) as a fimction of time are shown in Figs. 6, 7 and 8 for interictal, preictal and non- 
epileptic subjects, respectively. If SPR is close to the average under the null hypothesis (the 
broken line in the middle of the graph), then the null hypothesis cannot be rejected. 

30 However, if SPR, is too low, the null hypothesis must be rejected and the alternative 
hypothesis that ^ladj > ^remote (i e., Q2>0) must be accepted. Fig. 6 shows the SPR for a 
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collection of 40 interictal epochs of 20 minute duration taken from 8 patients with TLB. 
Each point represents one 40 second window. Fig. 7 is the SPR test for 24 one hour preictal 
epochs prior to seizure onset from this same cohort of patients, and Fig. 8 is the SPR for 24 
epochs of 20 minute duration taken from 6 non-epileptic subjects. 
5 Fig. 6 shows that the null hypothesis that the 52*s are the same for the adjacent and 

remote electrodes interictally must be rejected and the alternative hypothesis that Q2>0 must 
be accepted. Fig. 7 shows that up to about forty minutes prior to seizure onset, the 52's for 
the adjacent electrodes are significantly greater than those for the ipsilateral remote and the 
null hypothesis is below 5% significance level over that time period. Within about 40 
10 minutes prior to seizure onset, the SPR increases making rejection of the null hypothesis no 
longer possible. The marginal predictability of the electrodes adjacent to the site of ictal 
onset is significantly greater than that of the ipsilateral occipital electrodes, except within 
about half an hour prior to a seizure, at which time the marginal predictabilities take on 
similar values. 

15 Similar analyses comparing sphenoidal and occipital scalp electrodes have been 

preformed. The sphenoidal electrodes are also relatively close to the site of ictal onset in 
TLE. The results are qualitatively similar to those in Figs. 5 and 6, but are somewhat less 
distinct, probably due to increased noise and artifact from the sphenoidal electrodes. By 
comparison. Fig. 8 shows the summed-rank test for epochs derived from non-epileptic 

20 subjects. In marked contrast to the interictal epochs in patients with epilepsy. Fig. 8 shows 
that for non-epileptic subjects, the null hypothesis cannot be rejected, and the 52's for the 
ipsilateral temporal and occipital electrodes are statistically the same. 

The invention contemplates that the statistical results obtained using the present 
methods are not related to the behavior state of the epochs in the set. For example, test were 

25 run on the subsets of the preictal set and the interictal data for which the subjects were 

asleep during all the epochs, qualitatively, the predictive value of the present methods and 
systems was unchanged. 

Additional studies indicate that the results described in the previous paragraphs can 
be disaggregated and that the same effect can be seen in results from individual patients. 

30 For example. Fig. 9 shows show the values of Q2 for 4 interictal and 3 preictal epochs for 
one patient. It is clear that there is a systematic difference between the values of Q2 
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obtained for interictal as opposed to preictal epochs. In fact, a survey of the individual 
results for the eight patients studied reveals that, although for each given patient, Q2 during 
the preictal epoch is not necessarily statistically zero, it is always smaller than the typical 
values of Q2 during the interictal epochs for that patient. 
5 The patient results obtained from the methods and systems of the present invention 

strongly indicate a systematic change in a nonlinear measure computed on scalp EEG 
recordings prior to ictal onset in patients with epilepsy {e.g., medically refractory temporal 
lobe epilepsy). 

10 EXAMPLES 

The following examples are provided to demonstrate and further illustrate certain 
preferred embodiments of the present invention and are not to be construed as limiting the 
scope thereof. 

IS Example 1 

Subject Selection 

This example describes the selection processes used for selecting patients suitable 
for studies used to validate the methods of the present invention. Patients were evaluated by 
epileptologists at Henry Ford Hospital in Detroit. Presxu^gical evaluations followed a 

20 standardized protocol (A.M. Valachovic et aL, Language and its management in the surgical 
epilepsy patient in Medical Speech-Language Pathology, A.F. Johnson and B.H. Jacobson 
eds. Thieme, New York, NY, pp. 425-466 [1998]). In order to provide a homogenous group 
of patients, selection was limited to those afflicted with medically refractory mesiobasal 
temporal lobe epilepsy. Patients with this type of epilepsy are generally regarding as the 

25 most suitable candidates for epilepsy surgery. 

The specific criteria for inclusion in our analysis are: 1) seizures had to be of 
unilateral mesiobasal temporal lobe origin, documented by history, and interictal and ictal 
EEG recordings; 2) patients had to be between 18 and 60 years to reduce the likelihood of 
age related disorders such as cerebrovascular disease; 3) no mass lesions detectable with 

30 magnetic resonance imaging (MRI); 4) intelligence quotient of 70 or more; 5) no evidence 
of progressive neurological disorders, active neurological disorders other than epilepsy, and 
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no other significant medical disorder, severe depression or psychosis; 6) no evidence of 
damage to the hippocampus contralateral to the seizure focus as determined by MRI; 7) no 
history of substance abuse; 8) patients receiving barbiturates or benzodiazepines were 
excluded with the exception of intravenous benzodiazepines used for acute seizure control; 
5 and 9) no history of drug use other than antiepileptic drugs during the two weeks prior to the 
recordings. 

Example 2 
Electroencephalogram (EEG) Recordings 

1 0 Patient EEG recordings were recorded on a 128-channnel BMSI/Nicolet 5000 

System. (Nicolet Biomedical, Madison, WI). The band pass is 0.5 Hz to 100 Hz. The 
digital data is then transferred to a Linux workstation for conversion to ASCII text data and 
further analysis. An experienced epileptologist and a clinical neurophysiologist reviewed all 
EEG recordings. EEG recordings from the patients were visually inspected to identify 

15 epochs of interest for analysis. Epochs were divided into the following sets: 1) interictal, 
meaning at least 1 hour before and at least one hour after a seizure; 2) preictal, meaning 
within the hour preceding a seizure, and at least 1 hour following a seizure; and 3) ictal. 
Epochs were separated by behavioral state into: 1) wakefulness; 2) drowsiness; 3) stage 2 
non-REM sleep; 3) slow wave sleep; and 4) REM sleep. Waking and sleeping EEG from 

20 normal age and sex-matched subjects were analyzed. 

Example 3 

Study of MP and Different Behavior States 

In the following example an experienced epileptologist reviewed the complete scalp 
25 EEG (26 channels) for 61 interictal and 33 preictal epochs, each 20 minutes long from 14 
patients. The epileptologist categorized patient behavior during a plurality of 30 second 
interval of the epochs, and placed each interval into one of the following categories: 

Awake, eyes open — ^AEO 
30 Awake, eyes closed — ^AEC 
Lightly drowsy — ^Dl 
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Heavily drowsy — D2 

Stage 2 NonREM sleep — S2 

Stage 3 and 4 of NonREM sleep — S3/4 

REM sleep— REM 

5 

From the set of 40 thirty-second intervals for a given epoch, a summary behavior 
score for that epoch was produced. If 32 or more of the 30 second intervals (80%) of a 
given epoch were in the same behavior state, then that 20 minute epoch was deemed 
primarily in that behavior state (e.g., AEO or D2). If 60-79% of an epoch was spent in one 
10 state, then that epoch was considered as predominantly that state and indicated with the 
prefix P, thus PAEO or PD2. When less than 60% of an epoch was spent in just one state, 
then the epoch was considered a blend of 2 or more states, thus AEO/D2. Listed below are 
the numbers of various behavior states that comprise the basic data set of interictal and 
preictal epochs. 

15 

AEO— 49 
PAEO— 3 
D2— 9 
S2— 17 
20 Mixed states — 16 



No states were observed that were purely or predominantly Dl, AEC, S3/4 or REM, 
However, the states Dl and AEC do contribute to some of the mixed states. 

The data was used to test for the dependence of <Q2>, the value of Q2 averaged over 

25 20 minute epochs, on behavior state as well as on whether the epoch in question is preictal 
or interictal, and on whether the seizure focus is on the left or right side of the brain. The 
same methods were also used to test for the dependence of <62> and the marginal 
predictabilities for channels both remote from and adjacent to the seizure focus on the same 
set of variables. As the models become more complicated the more behavior state variables 

30 are included. For behavior states for which there are few observations, it is not advisable to 
include additional variables. There are two approaches to this issue. The first approach is 
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to agglomerate related behavior states (eg., AEO and AEC) into one category. The second 
approach is to eliminate those states with few observations. Two separate analyses were 
done using each of these approaches. In both cases, the number of behavior states against 
which the invention tested for dependence in either case was four (e.g., AEO, Dl, D2 and 
S2). For the first analysis, the invention used all 94 epochs of the basic data set. This set 
was called "inclusive.'' In this approach, the present invention categorized AEC 
observations with AEO observations. Note that the AEC observations only occurred in 
mixed states. The present invention also performed statistical tests on a subset of epochs, 
and their associated behavior states. However, with only a few representatives were 
eliminated. For this subset, all mixed states were removed. All epochs used in this data set 
were either purely or predominately one behavior state. This data set was called 
"restricted." The total number of epochs that comprise the restricted data was 78 and 
consisted of 54 interictal and 24 preictal epochs. The qualitative conclusions were the same 
for both data sets. 

To test for the dependence of <Q2> and < 5 2> on behavior state, on seizure focus 
location, and on preictal versus interictal, the present invention used a suite of Unear 
statistical models with up to five dummy variables. In this regard, the following (generally 
binary) variables were introduced: 

1) Xi=l if the behavior state is AEO and Xi=0 otherwise; 

2) X2=l if the behavior state is Dl, and X2=0 otherwise; 

3) X3=l if the behavior state is D2, and X3=0 otherwise; 

4) X4=l if the seizure focus is on the left side and X4=0 if the seizure focus is on the right 
side; and 

5) X5=l if the epoch is interictal and X5=0 if the epoch is preictal. 

Since the invention used four categories of behavior state, the assignment 
(Xi,X2,X3)=(0,0,0) uniquely corresponds to the behavior state S2. An exception to the 
assignment of binary values to the Xi, X2 and X3 occurred in the characterization of mixed 
states for the inclusive data set. In this case, the Xi took on values that reflected the firaction 
of the epoch associated with different behaviors. For example, a mixed state containing 
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25% of each of the four states, AEO, Dl , D2 and S2, was represented by the assignment 
X,=X2=X3=0.25. 

In the case of <Q2>, for example, the present invention constructed statistical tests 
for the dependence of <Q2> on the preictal versus interictal states, on the location of the 
5 seizure focus (i.e., left or right side of the subject's brain), and on the behavior state {see the 
following linear relation) 

(02 > = «o + Z ^J^J > (Fomiula 7) 

J 

where the aj are real coefficients and the Xj are a subset of the binary state variables defined 
10 above. By choosing different subsets of the binary variables, the present invention can test 
for the dependence of <Q2> on different behavioral, locational, and temporal states. For 
example. Formula 7 with a4=a5=0, can be used to test for the null hypothesis (/.e., that <Q2> 
is independent of behavior state). Similarly, ai=a2=a3=a5=0 can be set and the null 
hypothesis tested such that <Q2> does not depend on which hemisphere contains the seizure 
15 focus. The present invention can also set ai=a2=a3=a4=0, to test for the null hypothesis that 
<Q2>does not depend on whether the epoch is preictal or interictal. There is no indication 
that <Q2> depends on the side of the brain that contains the seizure focus, nor any indication 
that <Q2> depends on the behavior state. <Q2> does depend on whether the epoch in 
question is preictal or interictal. 
20 The present example also shows that analogous tests applied to the 20 minute 

averages of the MP's, < 82>, imply that there is no dependence on location of the seizure 
focus, nor on whether the epoch in question is preictal or interictal, but there is dependence 
on behavior state. 

Since <Q2> or < 8 2> do not necessarily have a linear dependence on state variables, 
25 the Unear fomi in Equation 7 was used as a basis to reject the null hypothesis that <Q2> (or 
< 62>) is independent of various state variables. It is possible that a more complicated form 
of dependence other than that of Equation 7 would reveal dependencies on state variables if 
the present approaches failed to reject the null hypothesis. One possibility, for example, is 
that a finer initial set of behavior states might reveal dependences. 

30 
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Results for Q2 

To perform the present analyses, the example provides a suite of six linear models of 
the form Equation 7. 

5 = ^0 + djATj (Test I) (Equation 8) 

An estimate of the value of as can be used to test for dependence of <Q2> on whether the 
epoch is interictal or preictal. The p-value for the estimate of as is 0.02 when the inclusive 
data set was used, and is 0,045 when restricted data set was used. P-values less than 0.05 

10 represent a rejection with greater than 95% confidence. Thus, the null hypothesis that a5=0 
can be rejected, which means that the null hypothesis that <Q2> was independent of whether 
the epoch in question was preictal or interictal can also be rejected. Moreover, the estimate 
of as was positive so that <Q2> for interictal epochs was significantly larger than it was for 
preictal epochs. This result was consistent with findings that show the dependence of <Q2> 

15 on the temporal proximity of a seizure. 

The dependence of <Q2> on seizure focus (e.g., left or right side of the brain) was 
tested. Two models were considered: ^ 

{Q2 ) = ^0 + ^4^4 (Test n) (Equation 9A) 

20 (Q2 ) = flo «4^4 + «5^5 (Test m) (Equation 9B) 



Statistical tests were constructed to reject the null hypothesis that 84=0. For the model 
Equation 9A, p-values of 0.22 and 0.24 for the inclusive and restricted data sets, 
respectively, were found. 

25 In Equation 9B the present invention used the partial F-test to test separately for the 

null hypotheses that as^O and that a4=0. It was found that the null hypothesis, 34=0 (p-value 
is 0.24 for the inclusive data set and 0.15 for the restricted data set), could not be rejected 
but that the null hypothesis, a5=0 (p-value is 0.03 1 for the inclusive data set and 0.032 for 
the restricted data set), could be rejected. These results are consistent with dependence of 

30 <Q2> on whether the epoch was interictal or preictal, but not dependent on seizure location. 
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The present invention also tested the dependence of <Q2> on behavior state, by 
considering two additional models: 

(Q2 ) = ao + a,^, + a^X^ -f a^X (Test IV) (Equation lOA) 

5 

(G2 ) = «o + ^5-^5 + + «2-^2 + «3^3 (Test V) (Equation lOB) 

In both of these models, the present invention tested the null hypothesis that ai=a2=a3=0. 

The alternative hypothesis was that at least one of the aj, j=l,2,3 was non-zero. For 
10 Equation lOA, the null hypothesis (the p-value is 0.11 for the inclusive data set and 0.13 for 

the restricted data set) could be rejected. To test the null in Equation lOB, the present 

invention used the partial F-test, again the null hypothesis could not be rejected (P-values 

are 0.1 8 for both the inclusive and restricted data sets). Consequently, <Q2> does not 

depend on behavior state. 
15 In another embodiment, the present invention provides a model that incorporates all 

five state variables simultaneously: 

5 

(G2> = «o +Z^y^j (Test VI) (Equation 11) 

20 The null hypothesis, ai=a2=a3=a4=0, was tested against the alternative hypothesis that at 

least one of the aj for j=l, 2,3,4 is non-zero. Using both the inclusive and restricted data sets, 
the partial F test shows that the null hypothesis, (the p-value is approximately 0.22 for each 
data set) cannot be rejected, and that <Q2> depends either on behavior state or on the 
location of seizure focus. The results of the tests for <Q2> are presented in Table L 

25 
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Table 1 



Summary of dependence tests for <Q2> 





<Q2> inclusive data set 


<Q2> restricted data set 




Null 

hypothesis 


p- value 


Null 

hypothesis 


p-value 


Test 
I 


Reject 


0.020 


Reject 


0.045 


Test 

n 


Accept 


0.22 


Accept 


0.24 


Test 

m 


(Accept, 
Reject) 


0.24, 0.031 
for (34, a5)=0 


(Accept, 
Reject) 


0.15, 0.032 
for (a4, a5)=0 


Test 
IV 


Accept 


0.11 


Accept 


0.13 


Test 
V 


Accept 


0.18 


Accept 


0.18 


Test 
VI 


Accept 


0.22 


Accept 


0.22 



5 Results for 82 

The results described above indicate that <Q2> was not sensitive to a subject's 
behavior state, but did depend on whether the epoch in question was preictal or interictal. 
Tests of the adjacent and remote channels using both inclusive and restricted data sets were 
conducted as follows: 

10 

(82) = +^5^5 (Test I) (Equation 12) 

Null hypothesis a5=0, alternate hypothesis as 0. 

(S2 ) = «o + «4-^4 (Test n) (Equation 13) 

1 5 Null hypothesis 34=0, alternate hypothesis 34 0. 
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(5 2 ) = «o + «4-*'4 + «5 ^5 (Test m) (Equation 14) 

Null hypotheses 34=0 or a5=0 using partial F-tests. 

(52) = ^o ^a^X^ -^a^X (TestIV) (Equation 15) 

Null hypothesis ai=a2=a3=0, alternate hypothesis, at least one of the aj (i=l,2,3,) is non-zero. 

(82) = ao + a^X^ + a\X^ + a^X^ + a^X^ (Test V) (Equation 16) 

Null hypothesis ai=a2=a3=0, alternate hypothesis, at least one of the a\ (i= 1,2,3,) is non-zero. 

5 

Test VI: (8 ^ ) = + ^ a.Xj (Test VI) (Equation 1 7) 

Null hypothesis ai=a2=a3=a4=0, alternate hypothesis, at least one of the aj (i=l. 2.3.4) is non- 
zero. 

The results for the tested of 82 are presented in the following tables. 

Table 2 



Test for dependencies of < d2> for a channel adjacent to the site of ictal onset 





<82> (adjacent) inclusive data set 


<62> (adjacent) restricted data set 




Null 

hypothesis 


p-value 


Null 

hypothesis 


p-value 


Test 
I 


Accept 


0.10 


Accept 


0.27 


Test 

n 


Accept 


0.47 


Accept 


0.58 


Test 

m 


(Accept, 
Accept) 


0.543. 0.50 for (a4, 
as)=0 


(Accept, 
Accept) 


0.50, 0.25 for (34, 
a5)=0 


Test 
IV 


Reject 


0.0003 


Reject 


0.0001 


Test 
V 


Reject 


0.00002 


Reject 


0.00002 


Test 
VI 


Reject 


0.00001 


Reject 


0.00007 
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Table 3 

Tests for dependencies of <62> for a channel remote from the site of ictal onset 





<S2> {remote) inclusive data set 


<82> (remote) restricted data set 




Null 

hypothesis 


p- value 


Null 

hypothesis 


p-value 


Test 
I 


Accept 


0.77 


Accept 


0.58 


Test 

n 


Accept 


0.62 


Accept 


0.64 


Test 

m 


(Accept, 

Accept) 


(0.77, 0.82) for (34, 

a5)=0 


(Accept, 
Accept) 


(0.60, 0.55) for 
(34, a5)=0 


Test 
IV 


Reject 


0.0001 


Reject 


0.0007 


Test 
V 


Reject 


0.00004 


Reject 


0.00003 


Test 
VI 


Reject 


0.00002 


Reject 


0.00009 



These results indicate that, unlike <Q2>, the individual <52> values from channels adjacent 
to and remote from the site of ictal onset strongly depend on behavior state. 



All publications and patents mentioned in the above specification are herein 
incorporated by reference. Although the invention has been described in connection with 
specific preferred embodiments, it should be understood that the invention as claimed 
should not be unduly limited to such specific embodiments. Indeed, various modifications 
of the described modes for carrying out the invention that are obvious to those skilled in the 
relevant fields are intended to be within the scope of the following claims. 
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